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1. Introduction and heuristics 



1.1 Evolution of spatially distributed systems in random media 

One of the often adequate and frequently used methods for studying the evolution of spatially dis- 
tributed systems under the influence of a random medium is homogenization. After rescaling, the 
system, modeled by partial differential equations with random coefficients, is approached by a sys- 
tem with 'properly averaged' deterministic coefficients, see e.g. |ZKU94j . But there are simple and 
important situations when random systems exhibit effects which cannot be recovered by such deter- 
ministic approximations and related fluctuation corrections. This concerns, in particular, localization 
effects for non-reversible random walks in random environment |S82j and for the electron transport in 
disordered media |And58j . 

Another such effect is that of intermittency. Roughly speaking, intermittency means that the so- 
lution of the system develops pronounced spatial structures on islands located far from each other 
that, in one or another sense, deliver the main output to the system. One of the sources of interest 
is magnetohydrodynamics and, in particular, the investigation of the induction equation with incom- 
pressible random velocity flelds |Z84j . jZMRS87j . Another source are simple mathematical models 
such as the random Fisher-Eigen equation that have been used to derive caricatures of Darwinian 
evolution principles |EEEF8^ . 

One of the simplest and most basic models exhibiting the effect of intermittency is the Cauchy 
problem for the spatially discrete heat equation with a random potential: 

dtu{t,x) = KAu{t,x) + £,{t,x)u{t,x), (t, j;) G (0, cx)) X Z'^, , , 

u{0,x) =uo{x), xel'^. ^ ' ' 

Here k > is a diffusion constant, A denotes the discrete Laplacian, 

y. |?;-x|=l 

^ is a space-time homogeneous ergodic random potential, and uq is a nonnegative initial function. 
Problem 1)1.1(1 is often called parabolic problem for the Anderson model or parabolic Anderson model 
(abbreviated PAM). As simplest localized initial datum one may take uq = 5q, and as non-local initial 
datum Uq = 1. In the latter case, the solution u{t, •) is spatially homogeneous and ergodic for each t. 
Let us remark that the solution u(t, x) of (|1.1|) allows the interpretation as average number of particles 
at site X at time t for branching random walks in random media given a realization of the medium ^, 
cf. |CM94j and the remarks in the next subsection. 

1.2 The PAM with time-independent potential 

In the particular case when the potential ^(t, x) = ^{x) is time- independent, the large-time behavior of 
the solution u to the PAM (|1.1() is determined by the spectral properties of the Anderson Hamiltonian 

n = KA + i (1.2) 

and therefore closely related to the mentioned localization of the electron transport. Namely, since 
(under natural assumptions on ^) the upper part of the spectrum of 7i in t^{U^) is a pure point 
spectrum |FMSS85j . |AM93j . the solut ion u admits the spectral representation 

u{t,-) = Y,(i^''{vi,uo)vi{-) (1.3) 

i 
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with respect to the random eigenvalues Aj and the corresponding exponentiahy locahzed random 
eigenfunctions Vi. (For simphcity we ignore the possible occurrence of a continuous central part of 
the spectrum.) As t increases unboundedly, only summands with larger and larger eigenvalues will 
contribute to (|1.3|) , and the corresponding eigenfunctions are expected to be localized more and more 
far from each other. Hence, for large t, the solution u{t, •) indeed looks like a weighted superposition 
of high peaks concentrated on distant islands. 

A mathematically rigorous understanding of the nature of the spectrum of the Anderson Hamiltonian 
and the random Schrodinger operator is still far from being complete. For an overview about some 
recent developments we refer to the surveys |BKSn4j and |LMWf)4| in this proceedings volume. The 
spectral results obtained so far do not yet seem directly applicable to answer the crucial questions 
about intermittency. A direct spectral approach clearly fails for space-time dependent potentials. 

In this survey we present a part of the results about intermittency for the PAM which have been 
obtained by use of more intrinsic probabilistic methods. In the next subsections we stick to the PAM 
with time-independent potential and localized initial datum: 

dtu{t,x) = KAu{t,x) + S,{x)u{t,x), {t,x) € (0,oo) x Z'^, , , 

u(0,x) = 5q{x), xeZ'^. ^ ' ^ 

We assume throughout that ^ = (^(a;))^^^^ is a field of i.i.d. random variables with finite positive 
exponential moments. Under these basic assumptions, u{t, x) has moments of all orders. 

The solution u to (|1.4|) describes a random particle flow in Z'^ in the presence of random sources 
(lattice sites x with ^(x) > 0) and random sinks (sites x with ^(x) < 0).^ Two competing effects are 
present: the diffusion mechanism governed by the Laplacian, and the local growth governed by the 
potential. The diffusion tends to make the random field u{t, •) flat, whereas the random potential ^ 
has a tendency to make it irregular. 

The solution u to H1.4|) also admits a branching particle dynamics interpretation. Imagine that 
initially, at time t = 0, there is a single particle at the origin, and all other sites are vacant. This 
particle moves according to a continuous-time symmetric random walk with generator kA. When 
present at site x, the particle is split into two particles with rate C+(x) and is killed with rate ^-(x), 
where ^+ = {£,+ {x))^^id and ^_ = {S^^{x))^^i^d are independent random i.i.d. fields {C-{x) may attain 
the value oo). Every particle continues from its birth site in the same way as the parent particle, and 
their movements are independent. Put ^(x) = 6.+ {x) — ^-(x). Then, given ^_ and the expected 
number of particles present at the site x at time t is equal to u{t, x). Here the expectation is taken over 
the particle motion and over the splitting resp. killing mechanism, but not over the random medium 

A very useful standard tool for the probabilistic investigation of ()1.4() is the well-known Feynman- 
Kac formula for the solution n, which (after time-reversal) reads 

u{t,x) =Eo expj^ i{X{s))ds^5^{X{t)) , (t,x) G [0, oo) x Z'^, (1.5) 

where (^(■s))se[o,oo) is continuous-time random walk on Z"^ with generator kA starting at x G Z'^ 
under E^,. 

Our main interest concerns the large-time behavior of the random field u{t,-). In particular, we 
consider the total mass, i.e., the random variable 

[/(t) = ^ n(t, x) = Eo expj / ^(X(s))ds|], t > 0. (1.6) 



X 



^Sites X with ^(a;) = — oo may be allowed and interpreted as ('hard') traps or obstacles, sites with ^(x) £ {—oo, 0) are 
sometimes called 'soft' traps. 
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Note that U{t) coincides with the value u{t, 0) of the solution u to the parabolic equation H1.4|) with 
initial datum ^(O, •) = ^ instead of u{0, •) = Jq- One should have in mind that, because of this, our 
considerations below also concern the large-time asymptotics of u. 
We ask the following questions: 

(1) What is the asymptotic behavior of U{t) as t — > oo? 

(2) Where does the main mass of n(t, •) stem from? What are the regions that contribute most to 
C/(i)? What are these regions determined by? How many of them are there and how far away 
are they from each other? 

(3) What do the typical shapes of the potential ^(•) and of the solution u{t, •) look like in these 
regions? 

We call the regions that contribute the overwhelming part to the total mass U(t) relevant islands or 
relevant regions. The notion of intermittency states that there does exist a small number of relevant 
islands which are far away from each other and carry asymptotically almost all the total mass U{t) of 
u{t,-). See Section [L3l for details. 

This effect may also be studied from the point of view of typical paths X{s), s S [0,t], giving 
the main contribution to the expectation in the Feynman-Kac formula H1.6|) . On the one hand, the 
random walker X should move quickly and as far as possible through the potential landscape to reach 
a region of exceptionally high potential and then stay there up to time t. This will make the integral 
in the exponent on the right of H1.6|) large. On the other hand, the probability to reach such a distant 
potential peak up to t may be rather small. Hence, the first order contribution to U{t) comes from 
paths that find a good compromise between the high potential values and the far distance. This 
contribution is given by the height of the peak. The second order contribution to U{t) is determined 
by the precise manner in which the optimal walker moves within the potential peak, and this depends 
on the geometric properties of the potential in that peak. 

It is part of our study to understand the effect of intermittency for the parabolic Anderson model in 
great detail. We distinguish between the so-called quenched setting, where we consider u{t, •) almost 
surely with respect to the medium ^, and the annealed one, where we average with respect to ^. It 
is clear that the quantitative details of the answers to the above questions strongly depend on the 
distribution of the field ^ (more precisely, on the upper tail of the distribution of the random variable 
1^(0)), and that different phenomena occur in the quenched and the annealed settings. 

It will turn out that there is a universal picture present in the asymptotics of the parabolic Anderson 
model. Inside the relevant islands, after appropriate vertical shifting and spatial rescaling, the potential 
^ will turn out to asymptotically approximate a universal, non-random shape, V, which is determined 
by a characteristic variational problem. The absolute height of the potential peaks and the diameter 
of the relevant islands are asymptotically determined by the upper tails of the random variable ^^(0), 
while the number of the islands and their locations are random. Furthermore, after multiplication 
with an appropriate factor and rescaling, also the solution u(t, •) approaches a universal shape on these 
islands, namely the principal eigenfunction of the Hamiltonian kA + V with V the above universal 
potential shape. Remarkably, there are only four universal classes of potential shapes for the PAM in 
(|1.4I) . see Section |1] for details. 

For a general discussion we refer to the monograph by Carmona and Molchanov |CM94j . the lectures 
by Molchanov |M94j . and also to the results by Sznitman about the important (spatially continuous) 
case of bounded from above Poisson-like potentials summarized in his monograph JS98|. A discussion 
from a physicist's and a chemist's point of view in the particular case of trapping problems (see 
also Section [2.21 below^. including a survey on related mathematical models and a collection of open 
problems, is provided in |HW94j . A general mathematical background for the PAM is provided in 
|C4M9nj . 
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1.3 Intermittency 

As before, let u denote the solution to the equation in H1.4|) with initial datum -^(0, •) = 1, but now 
with a homogeneous ergodic potential ^ = {^{x))^^^'^. Assume that all positive exponential moments 
of ^(0) are finite. Let Prob (•) and (•) denote probability and expectation w.r.t. ^. 

A first, rough, mathematical approach to intermittency consists in a comparision of the growth of 
subsequent moments of the ergodic field u{t, ■) as t ^ oo. Define 

Ap(t) = log(n(t,On, pGN, 

and write / < ^ if limt^oob(*) - /(*)] = 

Definition 1.1. For p G N \ {1}, the homogeneous ergodic field •) is called p-intermittent as 
t — > oo, if 

^ « ^. (1.7) 
p — 1 p 

Note that, by Holder's inequality, always Ap_i/(p — 1) < ^p/p- If the finite moment Lyapunov 
exponents 

Xp = lim -Ap(t), p£N, 

t — ^oo t 

exist, then the strict inequality Xp-i/{p — 1) < Xp/p implies p-intermittency. Such a comparison of 
the moment Lyapunov exponents has first been used in the physics literature to study intermittency, 
cf. [ZMRS87 , ZMRS88 . We will use this approach in Sectional 

To explain the meaning of Definition 11.11 assume (jl.Tf) for some p G N \ {1} and choose a level 
function ip such that Ap_i/(p — 1) <^ ip <^ ^p/p- Then, by Chebyshev's inequality, 

Prob('u(t,0) > e^^-W) < e-(P-i)^''W {u{t,Oy-'^) = exp{Ap_i(t) - {p - l)ip{t)} , 

and the expression on the right converges to zero as t — > oo. In other words, the density of the 
homogeneous point process 

r(t) = {x GZ"': u{t,x) > e^fW} 
vanishes asymptotically as t ^ oo. On the other hand, 

(u{t,Oyii^u{t,0) <e^''W}^ <eP^fW = gP^W-ApW {u{t,0)P) = o {{u{t,0)P)) 
and, consequently, 

(n(t,0)P) ~ (n(t,Of ]l{u(t,0) > e^^'W}) 
as i — > oo. Hence, by Birkhoff's ergodic theorem, for large t and large centered boxes B in Z'^, 

\B\''^Y,iL{t,xy \B\-^ u{t,xy. 

x€B x£Bnr{t) 

This means that the p-th moment {u{t, 0)^) is 'generated' by the high peaks of n(t, •) on the 'thin' set 
r(t) and therefore indicates the presence of intermittency in the above verbal sense. Unfortunately, 
this approach does not reflect the geometric structure of the set T(t). This set might consist of islands 
or, e.g., have a net-like structure. 

Theorem 1.2. //^ = (Cix)) xei,'' ^ non- deterministic field ofi.i.d. random variables with (e^^^^^) < 
oo for all t > 0, then the solution u{t, •) is p-intermittent for all p \ {1}. 

This is part of Theorem 3.2 in |GM90j . where, for general homogeneous ergodic potentials ^, neces- 
sary and sufficient conditions for p-intermittency of u{t, •) have been given in spectral terms of the 
Hamiltonian (|1.2|1 . 
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1.4 Annealed second order asymptotics 

Let us discuss, on a heuristic level, what the asymptotics of the moments of U{t) are determined by, 
and how they can be described. For simplicity we restrict ourselves to the first moment. 

The basic observation is that, as a consequence of the spectral representation ()1.3() . 

U{t) ^ e*^'(«) (1.8) 

(in the sense of logarithmic equivalence), where Xt{(p) denotes the principal (i.e., largest) eigenvalue 
of the operator kA + ip with zero boundary condition in the 'macrobox' Bt = [—t, t]'^ n Z'^. Hence, we 
have to understand the large-time behavior of the exponential moments of the principal eigenvalue of 
the Anderson Hamiltonian TC in a. large, time-dependent box. 

It turns out that the main contribution to {e''^^^^^) comes from realizations of the potential ^ having 
high peaks on distant islands of some radius of order a{t) that is much smaller than t. But this implies 
that At(^) is close to the principal eigenvalue of Ti on one of these islands. Therefore, since the number 
of subboxes of Bt of radius of order a{t) grows only polynomial in t and ^ is spatially homogeneous, 
we may expect that 



for R large as t ^ oo. 

The choice of the scale function a{t) depends on asymptotic 'stiffness' properties of the potential, 
more precisely of its tails at its essential supremum, and is determined by a large deviation principle, 
see ()1.15|) below. In Section|^we shall see examples of potentials such that a{t) tends to 0, to cxd, or 
stays bounded and bounded away from zero as t ^ oo. In the present heuristics, we shall assume that 
a{t) — > oo, which implies the necessity of a spatial rescaling. In particular, after rescaling, the main 
quantities and objects will be described in terms of the continuous counterparts of the discrete objects 
we started with, i.e., instead of the discrete Laplacian, the continuous Laplace operator appears etc. 
The following heuristics can also be read in the case where a{t) = 1 by keeping the discrete versions 
for the limiting objects. 

The optimal behavior of the field ^ in the 'microbox' is to approximate a certain (deter- 

ministic) shape (/? after appropriate spatial scaling and vertical shifting. It easily follows from the 
Feynman-Kac formula (|1.6() that 



H{t) 



where 



/?(*) = log(e*«W>, t>0, (1.9) 

denotes the cumulant generating function of ^(0) (often called logarithmic moment generating func- 
tion). Hence, the peaks of ^(•) mainly contributing to {U{t)) have height of order H{t)/t. Together 
with Brownian scaling this leads to the ansatz 

U)=c^itf[c{[-»m-^l (1-10) 

for the spatially rescaled and vertically shifted potential in the cube Qr = {—R, RY- Now the idea is 
that the main contribution to {U{t)) comes from fields that are shaped in such a way that ~ (/9 in 
Qij, for some ^p: Qr — > M, which has to be chosen optimally. Observe that 

It^^ in Qr ^ ^(.) ~ ^ + _^^(_^) ini?^,„(i). (1.11) 

Let us calculate the contribution to {U{t)) coming from such fields. Using (|1.8|) . we obtain 

{U{t) ^ in Qr}) « e^(*) exp{tA^„(t) {^^^i^)) } Prob ^^in Qr). (1.12) 
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The asymptotic scaling properties of the discrete Laplacian, A, imply that 

where A^((/?) denotes the principal eigenvalue of kA^ + lp in the cube Qr with zero boundary condition, 
and is the usual 'continuous' Laplacian. This leads to 

{U{t) ^ifin Qr}> ^ e^W e^p{-^XU^)} Prob ^ ip in Q^). (1.14) 

In order to achieve a balance between the second and the third factor on the right, it is necessary 
that the logarithmic decay rate of the considered probability is t/a{t)'^. One expects to have a large 
deviation principle for the shifted, rescaled field, which reads 

Prob {Ct ~ in Qr) exp|---^//j(v?)| , (1.15) 

where the scale a{t) has to be determined in such a way that the rate function Ir is non-degenerate. 
Now substitute (jLlSf) into 1)1.14(1 . Then the Laplace method tells us that the exponential asymptotics 
of {U{t)) is equal to the one of {U{t) I K, if m. Qr\ with optimal (p. Hence, optimizing on (p and 
remembering that R is large, we arrive at 

([/(t)>-e^Wexp{--^x}, (1.16) 

where the constant x is given in terms of the characteristic variational problem 

X=lim -uii [Ir{^)-X'r{^)]. (1.17) 

The first term on the right of (|1.16() is determined by the absolute height of the typical realizations 
of the potential and the second contains information about the shape of the potential close to its 
maximum in spectral terms of the Anderson Hamiltonian Ti in this region. More precisely, those 
realizations of ^ with ~ in Qr for large R and tp^ a minimizer in the variational formula in 
(|1.17() contribute most to {U{t)). In particular, the geometry of the relevant potential peaks is hidden 
via X ill the second asymptotic term of {U{t)). 



1.5 Quenched second order asymptotics 

Here we explain, again on a heuristic level, the almost sure asymptotics of U{t) as t ^ oo. Because of 
(|1.8|) . it suffices to study the asymptotics of the principal eigenvalue Xt{0- 

Like for the annealed asymptotics, the main contribution to Xt{£,) comes from islands whose radius is 
of a certain deterministic, time-depending order, which we denote a{t). As t — > oo, the scale function 
a{t) tends to zero, one, or oo, respectively, if the scale function a{t) for the moments tends to these 
respective values (see also ()1.2U|) below). However, a{t) is roughly of logarithmic order in a{t) if 
a{t) — > oo, hence it is much smaller than a{t). 

The relevant islands ('microboxes') have radius Ra{t), where R is chosen large. Let z € Bt denote 
the (certainly random) center of one of these islands B = z + Br^i^i) meeting the two requirements (1) 
the potential ^ is very large in B and (2) ^ has an optimal shape within B. This is further explained 
as follows. Let ht = maxst ^ be the maximal potential value in the large box Bt- (Then ht is a priori 
random, but well approximated by deterministic asymptotics, which can be deduced from asymptotics 
of H(t).) Then — ht is roughly of finite order within the relevant 'microbox' B. Furthermore, £, — ht 
should approximate a fixed deterministic shape in B. Hence, we consider the shifted and rescaled field 
in the box B, 

U-)=a{tfk{z + -a{t))-ht], in QR = {-R,Rf. (1.18) 
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Note that 

It^V inQR ^ ^(^ + .) + in^-z. (1.19) 

A crucial Borel-Cantelli argument shows that, for a given shape with probabihty one, for any 
t sufficiently large, there does exist at least one box B having radius Ra{t) such that the event 
{it ~ in Qr\ occurs if Ir{(p) < 1, where Ir is the rate function of the large deviation principle in 
H1.15() . If Ir{(p) > 1, then this happens with probability 0. For the Borel-Cantelli argument to work, 
one needs the scale function a(t) to be defined in terms of the annealed scale function a{t) in the 
following way: 

i.e., a{t) is the inverse of the map t i— > t/a{t)'^, evaluated at dlogt. Note that the growth of a{t) is 
roughly of logarithmic order of the growth of a{t), i.e., if the annealed relevant islands grow unbound- 
edly, then the quenched relevant islands also grow unboundedly, but with much smaller velocity. 

Hence, with probability one, for all large t, there is at least one box B in which the potential looks 
like the function on the right of (|1.19|) . The contribution to XtiO coming from one of the boxes B is 
equal to the associated principal eigenvalue 

^B-Aht + S(IF^(s(t))) ^ ^* + (l-^l) 

where we recall that Xr{^p) is the principal Dirichlet eigenvalue of the operator kA^ + 93 in the 'con- 
tinuous' cube Qr. Obviously, A((^) is asymptotically not smaller than the expression on the right of 
1)1.21(1 . In terms of the Feynman-Kac formula in (|1.5() . this lower estimate is obtained by inserting the 
indicator on the event that the random path moves quickly to the box B and stays all the time until 
t in that box. 

It is an important technical issue to show that, asymptotically as t — > cxd, A((^) is also estimated 
from above by the right hand side of (|1.21|) . if ip is optimally chosen, i.e., if A^((/9) is optimized over 
all admissible (p and on R. This implies that the almost sure asymptotics of U {t) are given as 

^ log U{t) « At(0 -h- t ^ 00, (1.22) 

t a(t)^ 

where x is given in terms of the characteristic variational problem 

X = lim inf [-X%{^)] . (1.23) 

R-*ooip: Qii^R,Iji(ip)<l 

This ends the heuristic derivation of the almost sure asymptotics of U{t). Like in the annealed case, 
there are two terms, which describe the absolute height of the potential in the 'macrobox' Bf, and the 
shape of the potential in the relevant 'microbox' B, more precisely the spectral properties of kA + S^ in 
that microbox. The interpretation is that, for R large and ip^, an approximate minimizer in (|1.23|) . the 
main contribution to U{t) comes from a small box B in Bt, with radius Ra{t), in which the shifted 
and rescaled potential looks like ip^. The condition lR{<^if) < 1 guarantees the existence of such a 
box, and A/j(99*) quantifies the contribution from that box. 

Let us remark that the variational formulas in (|1.23|) and ()1.17|) are in close connection to each other. 
In particular, it can be shown that the minimizers of (|1.23() are rescaled versions of the minimizers 
of (|1.17|) . This means that, up to rescaling, the optimal potential shapes in the annealed and in the 
quenched setting are identical. 



THE PARABOLIC ANDERSON MODEL 



9 



1.6 Geometric picture of intermittency 

In this section we explain the geometric picture of intermittency, still on a heuristic level. 

The heuristics for the total mass of •) in Section Fl. 51 makes use of only one of the relevant islands 
B in which the potential is optimally valued and shaped. In order to describe the entire function u{t, •), 
one has to take into account a certain (random) number of such islands. Let n{t) denote their number, 
and let zi, Z2, ■ ■ ■ , Zn(t) ^ Bt denote the centers of these relevant microboxes Bi, B2, ■ ■ ■ , whose 
radii are equal to Ra{t). Then, almost surely, 

n{t) 

U{t) = ^ u{t,x) u{t,x), ast^oo, (1.24) 

i.e., asymptotically the total mass of the random field u{t, •) stems only from the unions of the relevant 
islands, Bi, . . . , B^(j.y These islands are far away from each other. On each of them, the shifted and 
rescaled potential see (|1.18() ). looks approximately like a minimizer 99 of the variational problem 
in H1.23() . In particular, it has an asymptotically deterministic shape. This is the universality in the 
potential landscape: the height and the (appropriately rescaled) shape of the potential on the relevant 
islands are deterministic, but their location and number are random. 

The shape of the solution, u{t,-), on each of the relevant islands also approaches a universal de- 
terministic shape, namely a time-dependent multiple of the principal eigenfunction of the operator 
kA + (p^. 

2. Examples of potentials 



2.1 Double-exponential distributions 

Consider a distribution which lies in the vicinity of the double- exponential distribution, i.e., 

Prob(^(0) > r) exp{-e''/^}, r ^ cx), (2.1) 

with Q E (0, 00) a parameter. It turns out |GM98j that this class of potentials constitutes a critical 
class in the sense that the radius of the relevant islands stays finite as t ^ 00. This is related to the 
characteristic property of the double-exponential distribution that 

Prob(C(x) > /i) Prob(^(y) > h- Q\og2, ^{z) > h- Qlog2), 

meaning that single-site potential peaks of height h ^ 1 occur with the same frequency as two-site 
potential peaks with height of the same order. Hence, no spatial rescaling is necessary, and we put 
a(t) = 1. In Sections I3.1H,S.,SI below we shall describe our results for this type of potentials more 
closely. 

For the boundary cases g = 00 and g = ('beyond' and 'on this side of the double-exponential 
distribution, respectively), |(tM9 8I argued that the boundary cases a{t) J, and a(t) — > 00 occur. 
In other words, the fields beyond the double-exponential (which includes, e.g., Gaussian fields) are 
simple in the sense that the main contribution comes from islands consisting of single lattice sites. 
Unbounded fields that are in the vicinity of the case g = are called 'almost bounded' in [(irMQHj . 

2.2 Survival probabilities 

The case when the field ^ assumes the values —00 and only has a nice interpretation in terms 
of survival probabilities and is therefore of particular importance. The fundamental papers |DV75j 
and |DV79j by Donsker and Varadhan on the Wiener sausage contain apparently the first substantial 



10 



JURGEN GARTNER AND WOLFGANG KONIG 



annealed results on the asymptotics for the parabolic Anderson model. In the nineties, the thorough 
and deep work by Sznitman (see his monograph |S98j ). pushed the rigorous understanding of the 
quenched situation much further. 

2.2.1 Brownian motion in a Poisson field of traps. We consider the continuous case, i.e., the version of 
(|1.4|) with Z*^ replaced by and the lattice Laplacian replaced by the usual Laplace operator. The 
field ^ is given as follows. Let (xi)jg/ be the points of a homogeneous Poisson point process in W^, 
and consider the union O of the balls Ba{xi) of radius a around the Poisson points Xj. We define a 
random potential by putting 

= (2.2) 

^ ' \-oo if X G O. ^ ' 

(There are more general versions of this type of potentials, but for simplicity we keep with that.) The 
set O receives the meaning of the set of 'hard traps' or 'obstacles'. Let Tq = inf{t > 0: X{t) £ O} 
denote the entrance time into O for a Brownian motion (-'^(t))tg[o,oo)- Then we have the Feynman-Kac 
representation 

u{t, x) = Po {To > t, X{t) £ dx) jdx, 
i.e., u{t,x) is equal to the sub-probability density of X{t) on survival in the Poisson field of traps 
by time t for Brownian motion starting from the origin. The total mass U{t) = ¥q{To > t) is the 
survival probability by time t. It is easily seen that the first moment of U{t) coincides with a negative 
exponential moment of the volume of the Wiener sausage Use[o t] Ba{X{s)). 

Donsker and Varadhan analyzed the leading asymptotics of {U{t)) by using their large deviation 
principle for Brownian occupation time measures. The relevant islands have radius of order a{t) = 
^i/(d+2)_ To handle the quenched asymptotics of U{t), Sznitman developed a coarse-graining scheme 
for Dirichlet eigenvalues on random subsets of W^, the so-called method of enlargement of obstacles 
(MEO). The MEO replaces the eigenvalues in certain complicated subsets of by those in coarse- 
grained subsets belonging to a discrete class of much smaller combinatorial complexity such that 
control is kept on the relevant properties of the eigenvalue. 

Qualitatively, the considered model falls into the class of bounded from above fields introduced in 
Section 17!^ with 7 = 0. 

Related potentials critically rescaled with time have been studied in particular by van den Berg et 
al. fBSHnT] and by Merkl and Wiithrich ^MWVi2\ . 

2.2.2 Simple random walk among Bernoulli traps. This is the discrete version of Brownian motion 
among Poisson traps. Consider the i.i.d. field ^ = iS,{x)) ^^^.'^ where ^(x) takes the values or —00 
only. Again, u{t, x) is the survival probability of continuous-time random walk paths from to x 
among the set of traps O = {y £ Z'^ : (,{y) = —00}. 

In their paper |DV79j . Donsker and Varadhan also investigated the discrete case and described the 
logarithmic asymptotics of {U{t)) by proving and exploiting a large deviation principle for occupation 
times of random walks. Later Bolthausen |B94j carried out a deeper analysis of {U{t)) in the two- 
dimensional case using refined large deviation arguments. Antal |Ant94j . |Ant95j developed a discrete 
variant of the MEO and demonstrated its value by proving limit theorems for the survival probability 
U (t) and its moments. 



2.3 General fields bounded from above 

In |BK01aj and |BK01bj . a large class of potentials with esssup.^(0) < cxd is considered. Assume for 
simplicity that esssup ^(0) = and that the tail of ^ (0) at is given by 

Prob(C(0) > -x) PS expj-Dx^^l , x j 0, (2.3) 
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with D > and 7 € [0, 1) two parameters. The case 7 = contains simple random walk among 
Bernoulli traps as a particular case. The cumulant generating function is roughly H(t) ~ — const t"^, 
and the annealed scale function is a{t) ~ where v = (1 — 7)/(2 + d — ^7). The power v ranges from 
to l/{d + 2) as the parameter 7 ranges from 1 to 0. 

It turns out in |BK01aj that the rate function (see (|1.15() ) is given by 

lR{ip) = const / \ip{x)\~^ dx, (2.4) 
JQr 

where in the case 7 = we interpret the integral as the Lebesgue measure of the support of (p. The 
characteristic variational formula for the annealed field shapes in (jLlZf) has been analyzed in great 
detail in the case 7 = 0. In particular it was shown that the minimizer is unique and has compact 
support, and it was characterized in terms of Bessel functions. However, in the general case 7 G (0, 1), 
an analysis of (|1.17|) has not yet been carried out. 



2.4 Gaussian fields and Poisson shot noise 

Two important particular cases in the continuous version of the parabolic Anderson model are con- 
sidered in |GKOO) and ^GKMOO) (see also jCM95j for first rough results). The continuous version of 
()1.4() replaces Z"^ by and the discrete lattice Laplacian by the usual Laplace operator. Unlike in 
the discrete case, where any distribution on M may be used for the definition of an i.i.d. potential, in 
the continuous case it is not easy to find examples of fields that can be expressed in easily manage- 
able terms. Since a certain degree of regularity of the potential is required, the condition of spatial 
independence must be dropped. 

In |GK00j and |GKMOO] . two types of fields are considered: a Gaussian field ^ whose covariance 
function B has a parabolic shape around zero with B{0) = 0-2 > 0, and a so-called Poisson shot-noise 
field, which is defined as the superposition of copies of parabolic-shaped positive clouds around the 
points of a homogeneous Poisson point process in (in contrast to the trap case of Section 12.2(1 . 
A certain (mild) assumption on the decay of the covariance function (respectively of the cloud) at 
infinity ensures sufficient independence between regions that are far apart. 

Both fields easily develop very high peaks on small islands (the Poisson shot noise field is large 
where many Poisson points are close together). The annealed scale function is a{t) = t~^^'^ for the 
Gaussian field and a{t) = t'^/^e""''^'^ for the Poisson field |(;Knnj . 

3. Results for the double-exponential case 

In this section, we formulate our results on the large-time asymptotics of the parabolic Anderson 
model in the particularly important case of a double-exponentially distributed random potential, see 
Section 12.11 We handle the annealed asymptotics of the total mass U (t) in Section 13.11 the quenched 
ones in Section 13.21 and the geometric picture of intermittency in Section 13.31 The material of the 
first two subsections is taken from |(TM98j . that of the last subsection from ;GKMn4j . 



3.1 Annealed asymptotics 

As before, we assume that ^ = {£,{x))^^^d is a field of i.i.d. random variables. We impose the following 
assumption on the cumulant generating function of ^(0) defined by (|1.9() . 

Assumption (H). The function H{t) is finite for all t > 0. There exists g £ [0, cxd] such that 

^. H{ct)-cH{t) ^ , „ ,^ ^, 

lim = ^clogc for all a G (0, 1). 

t^oo t 
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Note that the vicinity of the double-exponential distribution H2.1|) corresponds to g £ (0, oo). If 
g = oo, then the upper tail of the distribution of ^(0) is heavier than in the double exponential case, 
whereas for ^) = it is thinner. 

Let 'P{Z'^) denote the space of probability measures on Z . We introduce the Donsker-Varadhan 
functional and the entropy functional on V{'L'^) by 

Sd{lj) = ^ {yji{x)- and Id{^J) = - ^ ^(x) log //(x). 



|a;-j;| = l 



respectively, and set 



Xd= inf + , £»g[0,oo]. (3.1) 

As before, let U{t) denote the total mass of the solution n(t, •) to the PAM (|1.4j) . 
Theorem 3.1. Let Assumption (H) he satisfied. Then, for any p G N, 

{U{tY) = exp {H{pt) - XdPt + o{t)} as t ^ oo. 

It turns out that Xd = "^dnxoiQ/K-), where xo- [0, oo) — > [0,1) is strictly increasing and concave, 
Xo(0) = 0, and Xo{q) ^ 1 as — > oo. Moreover, for g G (0, oo), each minimizer fj, of the variational 
problem (l^.lj) has the form fj, = const v'^ , where v = vi ^ ■ ■ ■ ^ Vd and each of the factor vi, . . . ,Vd is a 
positive solution of the equation 

kAv + 2gv log V = on Z 

with minimal ^^-norm. Uniqueness of v modulo shifts holds for large g/n but is open for small values 
of this quantity. 

Recall that U{t) = u{t,0), where u is the solution to (|1.4j) . but with non-localized initial datum 1 
instead of Jo- A much deeper question is the computation of the asymptotics of the 'correlation' 

_ {u{t,0)u{t,x)) 
(n(t,0)2) 

of the spatially homogeneous solution u of the PAM. Assuming additional regularity of the cumulant 
generating function H and uniqueness modulo spatial shifts of the minimizer of the variational 
problem 1)3. 1(1 . it was shown in ^GH9 9 that 

hm c{t, x) = — — . 

This indicates that the second moment (considered as lim/j^oo I-Bijl^-*^ Sx-s-Ba ""^(^j ^)) generated by 
rare high peaks of the solution u{t, •) with shape v{-). 



3.2 Quenched asymptotics 

Here we again consider i.i.d. random potentials {(,{x))^^^d in the vicinity of the double-exponential 
distribution H2.1() but formulate our assumptions in a different manner. 

To be precise, let F denote the distribution function of ^(0). Provided that F is continuous and 
F{r) < oo for all r € M (i.e., ^ is unbounded from above), we may introduce the non-decreasing 
function 

(^(r) = log ^ rGM. (3.2) 

1 — r (r) 
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Its left-continuous inverse -0 is given by 

i}{s) =mm{r £R: ip{r) > s}, s > 0. (3.3) 

Note that ip is strictly increasing with ip{^{s)) = s for all s > 0. The relevance of ip comes from 
the observation that has the same distribution as ^ o tj, where r] = {ri{x))^^^d is an i.i.d. field of 
exponentially distributed random variables with mean one. 

We now formulate our main assumption. 

Assumption (F). The distribution function F is continuous, F{r) < 1 for all r G M, and, in 
dimension d = 1, log \r\ F{dr) < oo. There exists g G (0, oo] such that 

liui[i{j{cs) - ip{s)] = glogc, cG(0,1). (3.4) 

If Q = OO, then ip satisfies in addition 

lim [V'(s + logs) - 0(s)] = 0. (3.5) 

s— >oo 



The crucial supposition 1)3. 4() specifies that the upper tail of the distribution of ^(0) is close to the 
double-exponential distribution ()2.1() for g G (0, oo) and is heavier for g = oo. Assumption ()3.5|) 
excludes too heavy tails. Note that 1)3. 5|) is fulfilled for Gaussian but not for exponential tails. 

The reader easily checks that 1)3. 4|) implies that 'ip{t) ~ glogt as t ^ oo. Let 

ht = max^(x), t > 0, (3.6) 

xeBt 

be the height of the potential ^ in St = [—t, tj'^nZ'^. It can be easily seen that, under Assumption (F), 
almost surely, 

ht = il^idlogt) + o{l) as t ^ oo. (3.7) 

Let us remark that it is condition H3.5|) which ensures that the almost sure asymptotics of ht in ()3.7|) 
is non-random up to order o(l). 

One of the main results in |GM98j . Theorem 2.2, is the second order asymptotics of the total mass 
U{t) defined in ([THIl . 

Theorem 3.2. Under Assumption (F), with probability one, 

log U{t) = t[ht-Xd + o(l)] as t ^ oo. (3.8) 

Here < Xd ^ 2dK. An analytic description of Xd is as follows. Define /: [— oo,0]^'' [0,oo] by 

„VO = ,('S^«'f^"". if. MO, CO), 

^ \{x £Z'^:V{x) > -oo}\, if^ = oo. ^ 



One should regard / as large deviation rate function for the fields i — ht (recall (|1.15j) and note that 
a{t) = 1 here). Indeed, if the distribution of ^ is exactly given by 1)2. then we have 



Prob -h> V{-) in Z'^) = exp |-e'^/^/(y)} 



for any V : Z"^ ^ [-oo,0] and any h G (0,oo). For V G [-oo,0]^ , let XiV) G [-oo,0] be the top 
of the spectrum of the self-adjoint operator kA -|- y in the domain {V > —oo} with zero boundary 
condition. In terms of the Rayleigh-Ritz formula, 

X{V)= sup {{KA + V)fj), (3.10) 

/e^2(Zd): ||/||2 = 1 
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where (•, •) and || • II2 denote the inner product and the norm in ^^(Z'^), respectively. Then 

-Xd = snp{x{V):V e[-^,0f, I{V) < l} . (3.11) 
This variational problem is 'dual' to the variational problem ()3.1() and, in particular, Xd = Xd- 

3.3 Geometry of intermittency 

In this section we give a precise formulation of the geometric picture of intermittency, which was 
heuristically explained in Section [L6l 

We keep the assumptions of the last subsection. For our deeper investigations, in addition to 
Assumption (F), we introduce an assumption about the optimal potential shape. 

Assumption (M). Up to spatial shifts, the variational problem in (j3.11|) possesses a unique maxi- 
mizer, which has a unique maximum. 

By we denote the unique maximizer of ()3.11|) which attains its unique maximum at the origin. 
We will call Ki, optimal potential shape. Assumption (M) is satisfied at least for large g/ k. This 
fact as well as further important properties of the variational problem 1)3. 11() are stated in the next 
proposition. 

Proposition 3.3. (a) For any g € (0, cxd], the supremum in (|3.11j) is attained. 

(b) // g/n is sufficiently large, then the maximizer in (|3.1H) is unique modulo shifts and has a 
unique maximum. 

(c) If Assumption (M) is satisfied, then the optimal potential shape has the following properties. 

(i) If Q ^ (Oj 00), then K = /* Cg) • • • ^ /* for some : Z — > (—00, 0). If g = 00, then K is 
degenerate in the sense that V*(0) = and V^{x) = —00 for x ^ 0. 

(ii) The operator kA + V^, has a unique nonnegative eigenfunction u;* G ^^(Z'^) with ii;*(0) = 1 
corresponding to the eigenvalue X{V^,). Moreover, G £^(Z^). If g & (0,oo), then is 
positive on 11^ , while = 5q for g = cc. 

We shall see that the main contribution to the total mass U{t) comes from a neighborhood of the set 
of best local coincidences of ^ — with spatial shifts of V^. These neighborhoods are widely separated 
from each other and hence not numerous. We may restrict ourselves further to those neighborhoods 
in which, in addition, n(t, •), properly normalized, is close to w^:. 

Denote by Bji{y) = y + Br the closed box of radius R centered at y G Z*^ and write 

Br{A) = J BR{y) (3.12) 

for the 'i2-box neighborhood' of a set A C Z'^. In particular, Bq{A) = A. 

For any e > and any sufficiently large g G (0, 00], let r(e, g) denote the smallest r G No such that 

11^*112 ^ w^{x) < e. (3.13) 

Note that r(e, cxd) = 0, due to the degeneracy of Woo. Given /: Z*^ — > R and > 0, let ||/||_r = 
sup^'eSfl \f{x)\- 

The main result of |(TKMn4! is the following. 
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Theorem 3.4. Let the Assumptions (F) and (M) he satisfied. Then there exists a random time- 
dependent subset r* = r*j^g2j '^f ^tiog^t such that, almost surely, 

(i) liminf-i- V u{t,x)>l-£, eG(0,l); (3.14) 

t^oo U[t) ^ — ' 

^eBr(.,e)(r*) 

(m) |r*|<t°(^) and ^min Jy - y\ > t^""^^'^ as t ^ oo; (3.15) 

(Hi) lim max||^(y + •) - /it - V;(-)L = 0, i? > 0; (3.16) 

{iv) limmax - «;,(•) =0, > 0. (3.17) 



t— >oo j/gr* 



u{t,y) 



Theorem 13.41 states that, up to an arbitrarily small relative error e, the islands with centers in 
r* and radius r(e, g) carry the whole mass of the solution u{t, •). Locally, in an arbitrarily fixed R- 
neighborhood of each of these centers, the shapes of the potential and the normalized solution resemble 
ht + V* and w*, respectively. The number of these islands increases at most as an arbitrarily small 
power of t and their distance increases almost like t. Note that, for g = oo, the set i?j,(j^g)(r*) in (|3.14l) 
is equal to F* and, hence, the islands consist of single lattice sites. 

It is an open problem under what assumptions on the potential ^ the number |r*| of relevant peaks 
stays bounded as t — > oo; we have made no attempt to choose T* as small as possible. 

4. Universality 

In this section we explain that, under some mild regularity assumptions on the tails of ^(0) at its 
essential supremum, there are only four universality classes of asymptotic behaviors of the parabolic 
Anderson model. Three of them have already been analyzed in the literature: the double-exponential 
distribution with g G (0, oo) respectively = oo GM98^ 



GH99 


. ,GKon,. 


tGKMOO 




GKM04 


) (see 


BKOTal. IBKnibl. RMI. 


IAnt94l. IAnt95l 



(see Section [2. 3|) . A fourth and new universality class is currently under investigation, see |IIKMn4] . 
This class lies in the union of the boundary cases = of the double-exponential distribution ('almost 
bounded' fields) and 7 = 1 for the general bounded from above fields. Examples of distributions that 
fall into this class look somewhat odd, but it turns out that the optimal potential shape and the 
optimal shape of the solution are perfectly parabolic, respectively Gaussian, which makes this class 
rather appealing. In particular, the appearing variational formulas can be easily solved explicitly and 
uniquely. 

We now summarize |HKM04] . Our basic assumption on the logarithmic moment generating function 
H in (|1.9() is the following. 

Assumption (H): There are a function H : (0, 00) M and a continuous auxiliary function 
rj: (0,00) (0,00) such that 

lim^M_^=#(,)^0 for, /I. (4.1) 



The function H extracts the asymptotic scaling properties of the cumulant generating function H. 
In the language of the theory of regular functions, the assumption is that the logarithmic moment 
generating function H is in the de Haan class, which does not leave many possibilities for H: 

Proposition 4.1. Suppose that Assumption (H) holds. 
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(i) There is a 7 > such that limt|oo v{yt)/'n{'^) — /o*" (^it-V y > 0, i.e., ij is regularly varying of 
index 7. In particular, ri{t) = t'>'+°(i) as t ^ 00. 



(ii) There exists a parameter p > such that, for every y > 0, 

(a) %)=/,^^ j/7/1, 

(b) H{y) = pylogy i/7 = 1. 

Our second regularity assumption is a mild supposition on the auxiliary function rj. This assumption 
is necessary only in the case 7 = 1 (which will turn out to be the critical case). 

Assumption (K): The limit r/* = lim^^oo ^(0/^ ^ [0, co] exists. 
We now introduce a scale function a: (0, cxd) — > (0, cxd), by 

tait)-d a{tY ■ ^ ' ' 

The function a{t) turns out to be the annealed scale function for the radius of the relevant islands in 
the parabolic Anderson model. We can easily say something about the asymptotics of a{t): 

Lemma 4.2. Suppose that Assumptions (H) and (K) hold. If j < I and r/^, < cxd, then there exists a 
unique solution a: (0, 00) (0,oo) to ()4.2() . and it satisfies limt^oo ta{t)~'^ = cxd. Moreover, 

(i) 7/7 = 1 and < ry* < cxd, then lim^^oo 0({t) = iZ-^/r/T G (0, 00). 

(ii) 7/7 = 1 and r]^ = 0, then a{t) = as t^ 00, where v = {I - -f) / {d + 2 - d-y). 

Now, under Assumptions (H) and (K), we can formulate a complete distinction of the PAM into 
four cases: 

(1) ?7* = 00 (in particular, 7 > 1). 

This is the boundary case ^> = 00 of the double-exponential case. We have a{t) ^ as t — > 00, 
as is seen from 1)4. 2() . i.e., the relevant islands consist of single lattice sites. 

(2) ?7* G (0,00) (in particular, 7 = 1). 

This is the case of the double-exponential distribution in Section 12.11 By rescaling, one can 
achieve that r/^, = 1. The parameter g of Proposition I4.1f ii')(b) is identical to the one in 
Assumption (H) of Section 13.11 

(3) ?7* = and 7 = 1. 

This is the case of islands of slowly growing size, i.e., a{t) — > cxd as t ^ 00 slower than any 
power of t. This case comprises 'almost bounded' and bounded from above potentials. This 
class is the subject of |HKM04j , see also below. 

(4) 7 < 1 (in particular, r]^ = 0) 

This is the case of islands of rapidly growing size, i.e., a{t) — > cxd as t ^ 00 at least as fast 
as some power of t. Here the potential ^ is necessarily bounded from above. This case was 
treated in |BK01aj : see Section [T^ 

Let us comment on the class (3), which appears to be new in the literature and is under investigation 
in |HKM04] . One obtains examples of potentials (unbounded from above) that fall into this class by 
replacing g in 1)2. 1|) by a sufficiently regular function g{r) that tends to as r ^ cxd, and other examples 
(bounded from above) by replacing 7 in 1)2. 3(1 by a sufficiently regular function 7(x) tending to 1 as 
X I 0. According to Lemma 14.21 the scale function a{t) defined in ()4.2|) tends to infinity, but is slowly 
varying. The annealed rate function for the rescaled potential shape, Ir, introduced in (jl.lSp turns 
out to be 

Ir{^) = const [ e-^^^^/^dx. (4.3) 
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The characteristic variational problem for the annealed potential shape in ()1.17|) . turns out to be 
uniquely minimized by a parabolic function 9J*(x) = const — ^3||x||2, and the principal eigenfunction 

II 1 1 2 

of the operator kA + (p^, is the Gaussian density = const e~^ll^ll2 . 

5. TiME-DEPENDENT RANDOM POTENTIALS 

In this section we study the intermittent behavior of the parabolic Anderson model (PAM) with a 
space-time homogeneous ergodic random potential : 

dt u{t, x) = KAu{t, x) + [^(t, x) - x))] u{t, x), (t, x) £ (0, oo) x Z"^, , , 

u{0,x) = 1, x£ Z'^. ^ ^ 

Note that for time-dependent potentials the direct connection to the spectral representation of the 
Anderson Hamiltonian (|1.2|) is lost. Our focus will be on the situation when the potential ^ is given 
by a field {Yk{t);k € N} of independent random walks on Z'^ with diffusion constant g in Poisson 
equilibrium with density v. 

e(t,x)=7E'^nw(^)' (5-2) 

fcGN 

where 7 denotes a positive coupling constant. Clearly 

{C{t,x))=uj. 

This deterministic correction to the potential in (|5.1() has been added for convenience to eliminate 
non-random terms. 

The form ()5.2() of the potential is motivated by the following particle model. Consider a system of two 
types of independent particles, A and B, performing independent continuous-time simple random walks 
with diffusion constants k and g and Poisson initial distribution with densities v and 1, respectively. 
Assume that the 5-particles split into two at a rate that is 7 times the number of j4-particles present 
at the same location and die at rate 1^7. Hence, the A- and i?-particles may be regarded as catalysts 
and reactants in a simple catalytic reaction model. Then the (spatially homogeneous and ergodic) 
solution u{t, x) of the PAM 1)5. 1|) is nothing but the average number of reactants at site x at time t 
given a realization of the catalytic dynamics (|5.2|) . Such a particle model (with arbitrary death rate) 
has been considered by Kesten and Sidoravicius _KS03j . We will come back to the results in |KS03j at 
the end of this section. We further refer to the overview papers by Dawson and Fleischmann |DF00j 
and by Klenke |K100j for continuum models with singular catalysts in a measure- valued context where 
questions different from ours have been addressed. 

Our aim is to study the moment Lyapunov exponents 

Xp= lim ilog(n(t,0)P) 

t— >oo t 

as well as the quantities 

A;= hm iloglog(n(t,0)^') 

{p = 1,2,...) as functions of the model parameters. The phase diagram will turn out to be different 
in dimensions d = 1,2, d = 3, and d > 4. 

Definition 5.1. a) For p G N, we will say that the PAM (|5.1|) - H5.2() is strongly p-catalytic if A* > 0. 
Otherwise the PAM will be called weakly p-catalytic. 

b) For p G N \ {1} and Ap < 00, we will say that the PAM is p-intermittent if Xp/p > Xp-i/{p — 1). 
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We believe that strongly catalytic behavior is related to heavy tails of the Poisson distribution of 
the catalytic point process {Yk{t); /c S N} and may occur if the main contribution to the p-th moment 
comes from realizations with a huge number of catalysts at the same lattice site. Recall that p- 
intermittency means that, for large t, the p-th moment, considered as limi^^oo SxgBh ^(^'^)^' 

is 'generated' by high peaks of the solution u{t, ■) located far from each other. 

For potentials of the form 

C{t,x)=jW,{t) 

with iWx{t)),j.^^d being a field of i.i.d. (or correlated) Brownian motions and (|5.1() understood as a sys- 
tem of ltd equations, the moment Lyapunov exponents have been shown by Carmona and Molchanov 
|CM94j to exhibit the following behavior. In dimensions d = 1,2 there is p-intermittency for all 
p G N \ {1} and all choices of the model parameters k and 7. If d > 3, then there exist critical 
points < C2 < C3 < • • • such that p-intermittency holds if and only if K/7 < Cp. For this model, the 
asymptotics of the almost sure ('quenched') Lyapunov exponent as k — > has been investigated in 
\rMM\ . |(:MV96j . ICKMOH . and |(:MSn2j . 

In the following we present the results for catalytic potentials of the form (|5.2)) obtained in |(THn4j . 
Our analysis of the moment Lyapunov exponents is based on the following probabilistic represen- 
tation of the p-th moment which is easily derived from the Feynman-Kac formula for u(t,0): 

(n(t, or) = eJ^) exp i^i^lj^Y. ^(*' ^^(*)) ' (5-3) 

where k'^^ denotes expectation with respect to p independent random walks Xi, . . . ,Xp on Z"^ with 
generator kA starting at the origin, and w is the solution of the random initial value problem 

dtw{t,x) = QAwit,x)+-f^P^^6x,(t){x)il + w{t,x)), {t,x) £ {0,00) xZ^, 

w{0,x) =0, x£ Z*^. 

One of the main difficulties in analyzing (|5.3|) is related to the observation that w{t, Xi{t)) depends 
in a nontrivial way on the whole past {Xj(s); < s < t}, j = 1, . . . ,p, of our random walks (although 
our notation does not reflect this). 

For r > 0, let //(r) denote the upper boundary of the spectrum of the operator A -|- r6o in £^(Z°'). 
It is well-known that, in dimensions d = 1,2, ^{r) > for all r and, in dimensions d >3, fi{r) = for 
< r < and /i(r) > for r > r^, where 

ra = 1/Gd{0) 

and Gd denotes the Green function associated with the discrete Laplacian A. 

Theorem 5.2. For any choice of the parameters of the PAM it5. the limit X* exists, and 

>^l = qKpi/q), pen. 

Hence, for d = 1,2, the PAM (|5.1() - (|5.2() is always strongly p-catalytic, whereas for d > 3 this is 
true only ii p'j/g exceeds the critical threshold r^,. Note that A* does not depend on k nor on u. 

We next study the behavior of the moment Lyapunov exponents Xp = Ap(k) as a function of the 
diffusion constant k in the weakly catalytic regime < p^/ g < r^ for d>3. We will mainly describe 
their behavior for small and for large values of k. 

Theorem 5.3. Let d > 3, p £ N, and < p^/g < r^. Then the limit Xp exists and is finite for all k 
and V. Moreover, k i— > Ap(«:) is strictly decreasing and convex on [0, 00) and satisfies 

Xp{K) Ap(0) pj/g 
lim — = — = — . 

Kio p p rd- pj/ g 
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Hence, in any dimension d> 3, the PAM ()5.1|) - H5.2|) is p-intermittent in the weakly catalytic regime 
for small values of the diffusion constant k. 

To formulate the behavior of the moment Lyapunov exponents for k — > oo, we introduce the 
variational expression V for the polaron problem analyzed in |L77j . |DV83j . and |BDS93j : 

2 



V 



sup 

fSC°° (k3 
Il/ll2 = l 



l|V/||^ 



Theorem 5.4. Let d > 3, p £ N, and < pj/g < r^. Then 

P i/ii>4. 

In other words, in dimensions d = 3 for p in the weakly catalytic regime, the PAM H5.1() ~ (|5.2|) is 
p-intermittent also for large values of n. We conjecture intermittent behavior for all k. In dimensions 
d > 4, the leading term in the asymptotics of \p{ti)/p as k — > oo is the same for all p in the weakly 
catalytic regime. We in fact conjecture that there is even no intermittency in high dimensions for 
large k. 

Let us remark that in |KS03j Kesten and Sidoravicius obtained the following results for the above 
mentioned catalytic particle model with arbitrary death rate 5 instead of v^. If d = 1,2, then for 
all choices of the parameters, the average number of -B-particles per site tends to infinity faster than 
exponential. If d > 3, 7 small enough and 5 large enough, then the average number of i?-particles per 
site tends to zero exponentially fast. The first of these two results is covered by Theorem 15.21 since 
the death rate 5 does not affect A^. The second result is covered by Theorem 15. 3| since < Ai < 00 in 
the weakly catalytic regime and the death rate 5 shifts Ai by — 5. Kesten and Sidoravicius further 
show that in all dimensions for 7 large enough, conditioned on the evolution of the >l-particles, there 
is a phase transition. Namely, for small 6 the i?-particles locally survive, while for large 5 they become 
locally extinct. In j(THn4j there are no results for the quenched situation. The analysis in |KSn3j 
does not lead to an identification of Lyapunov exponents, but is more robust under an adaption of 
the model than the above analysis based on the Feynman-Kac representation. 
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